Role of Inelastic Tunneling through the Barrier in Scanning Tunneling Microscope 

Experiments on Cuprates 
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The tunneling patli between tlie Cu02-layers in cuprate superconductors and a scanning tunneling 
microscope tip passes through a barrier made from other oxide layers. This opens up the possibility 
that inelastic processes in the barrier contribute to the tunneling spectra. Such processes cause one 
or possibly more peaks in the second derivative current-voltage spectra displaced by phonon energies 
from the density of states singularity associated with superconductivity. Calculations of inelastic 
processes generated by apical O-phonons show good qualitative agreement with recent experiments 
reported by Lee et al. U). Further tests to discriminate between these inelastic processes and coupling 
to planar phonons are proposed. 

PACS numbers: 74.72.-h,74.25.Kc,74.50.+r,68.37.Ef 



The importance of the clectron-phonon interaction in 
the high temperature cuprate superconductors has long 
been a matter of debate. In conventional superconduc- 
tors tunneling experiments analyzed by McMillan and 
Rowell not only unequivocally established that the 
exchange of phonons between electrons is the underlying 
mechanism for Cooper pairing but also allowed its spec- 
troscopic determination. In the case of cuprates there 
are many good reasons to doubt the dominance of an at- 
traction due to phonon exchange. These reasons range 
from the d-wave rather than s-wave pairing symmetry, 
to their electronic structure as lightly doped Mott insu- 
lators. The latter points toward a strong onsite Coulomb 
repulsion which cannot be treated perturbatively, vio- 
lating a key postulate in the standard BCS-Eliashberg 
theory. A number of proposals have been made for a 
subdominant role for the electron-phonon interactions 
particularly involving the half-buckling phonon modes 
which couple attractively in the d-wave pairing chan- 
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Very recently, a new set of tunneling measurements at 
low temperatures has been reported by Lee et al. Q us- 
ing a STM (Scanning Tunneling Microscope). Results 
with this local tunneling probe show considerable vari- 
ations in the form and energy of the density of states 
(DOS) singularity associated with the superconductivity 
and also a peak in the second derivative current-voltage 
spectra which is displaced from the DOS peak by a fairly 
constant energy. Lee et al. [J report further that iso- 
tope substitution on the 0-ions shifts this energy con- 
sistent with a phonon origin for this peak. This raises 
the question whether these spectra should be interpreted 
analogously to those of conventional superconductors or 
if these peaks have another source. 

A series of works explain the appearence of side peaks 
in the DOS of d-wave superconductors in terms of lo- 
cahzed bosonic modes , resonant spin modes 0, Q| 
or the breathing mode phonons 0|. All these propos- 
als consider inelastic processes that happen directly in 



the superconducting plane. In this letter we examine 
an alternative explanation in terms of inelastic tunnel- 
ing processes. Such inelastic processes are well estab- 
lished in electron tunneling spectroscopy [isl Iisl 
of metal-insulator-metal junctions and in STM spectra 
taken on metals with adsorbed molecules m [ig. In 
cleaved BSSCO the topmost layer is the BiO-plane, but 
the electrons tunnel into and out of the Cu02-planes. In- 
elastic tunneling can occur in the intervening barrier, e.g. 
when the tunneling hole passes through the apical 0-ion. 
Although the 2p2-orbitals of the apical 0-ion do not hy- 
bridize with the 3d2;2-j,2-orbitals of the Cu-ion directly 
below there is a weaker hybridization with orbitals cen- 
tered on neighboring Cu-ions However we will not 
go into such details here. Rather, we examine a simple 
model which has an important tunneling path through 
an apical 0-orbital removed from the chemical potential. 
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FIG. 1: Left panel: Geometry relevant for the scanning 
tunneling microscope experiment, the superconducting Cu02- 
plane lies below BiO- and SrO-layers. Right panel: Sug- 
gested inelastic cotunneling process via the apical oxygen 
atom leading to phonon-satellites in the current-voltage char- 
acteristics in electron-energy scheme. 



This tunneling path is in turn coupled to a vibration of 
this ion. Our aim is to establish the key features of this 
inelastic tunneling process and to compare them to the 
experimental spectra. 

Our aim is to calculate the cotunneling current through 
the 2p2-orbital of the apical 0-ion. For this purpose we 
describe the setup depicted in Fig.^by the Hamiltonian, 

H = H„ + iJs + Ht 
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where represents the normal conducting STM tip, 
Hg the superconducting Cu02-planc and Ht the hop- 
ping onto the apical oxygen orbital. Operator pi cre- 
ates a hole in the oxygen p-orbital, operator b^^ creates 
a phonon with momentum q. The STM tip is charac- 
terized by a constant DOS Nyi{e) oc nn(e) = 1. Partial 
line widths FjuFs of the 2pz-orbital at energy Eq are 
used to describe its hybridization with the STM tip and 
the superconducting Cu02-plane. As we shall discuss at 
the end of this section the above parameters mainly af- 
fect the overall prefactor in the tunneling current and 
are not very crucial. The shape of the current-voltage 
characteristics is determined by the DOS of the super- 
conductor and the phonon spectrum. For the super- 
conductor we assume for simplicity a d-wave type DOS 
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However, in 



N,{e) (xn^ie) = (Re[e/Ve 

the experiment [l| the shape of the superconducting DOS 
varies strongly from point to point. Since the phonon 
involves mainly the motion of the apical oxygen atoms 
we assume an optical phonon band centered around fre- 
quency r^o with a weak dispersion (i.e. a small band 
width W). 

The derivation of the cotunneling current follows 
closely Ref. which deals with a similar non- 

equilibrium situation in the normal conducting state. 
Tunneling is treated to lowest order, while electron- 
phonon coupling is treated non-perturbatively allowing 
for multi-phonon processes. Applying a Lang-Firsov 
transformation |21l and standard manipulations of many- 
body physics |23| one arrives at the following expression 
for the cotunneling current (we set h = I) 



2tt 
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The factors F.j;{e), Fx{e) with a: = n, s are obtained from 
a convolution, 



dujnx{£ + w) 
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B{-Lu)[l- Ue + u)] J ' 
(3) 

of the relevant DOS, nx{s), and the Fermi factors, fx{£), 
with the phonon correlation function B{uj). The lat- 
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FIG. 2: Shape of dl /dV for different widths W of the phonon 
band. Parameters are A = 38meV, fio = 50meV,A = 0.4, 
Eo = 2eV. 



ter will be discussed in detail below. The retarded self 
energy S](e) of the 2pz-orbital reads 

i 
2 



(4) 



If the energy Eq of the virtual oxygen state is the largest 
energy scale in the problem, the energy de pen dence of 
the denominator of Eq. ^ can be neglected [23j . In this 
case, we may approximate the tunneling current by 



I{V) ~ / {F^{e)F,{e) - F,{e)F^{e)}de 



(5) 



and the parameters i?o, Fn, Fg describing the virtual state 
enter only in the (unimportant) prefactor. Even if terms 
with e ~ Eq contribute significantly to the integral ||5J), 
the energy dependence of the denominator will only lead 
to an asymmetry of the current-voltage signal on a larger 
energy scale. In the following, we therefore concentrate 
on the influence of the phonon spectrum which creates 
the sharp features visible in the experiment P|. 

The phonon spectrum enters the current-voltage spec- 
trum (0) via the modified occupation factors Q. In 
Eq. 0, the phonon correlation fimction B{u!) = 
{2tt)-^ J dte"^'^{X{t)X'^{0)) describes possible phonon 
emission or absorption of the oxygen 2p2-orbital. Those 
processes lead to an effective shift of the Fermi occupa- 
tion factors of both the superconducting plane and the 
normal conducting STM tip. The "displacement" opera- 
tor X ~ cxp Mq(6g — 6Lg)/ijJg| results from a Lang- 
Firsov transformation |2 j| . Replacing the g-summation 
by a frequency integration over a box density of states 
and assuming zero temperature for the phonon bath we 
obtain the correlation function 



(X(i)Xt(0))=cxp^A^(^sin-^, 
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FIG. 3: Shape of dl/dV for different electron-pfionon cou- 
plings A. Parameters are A = 38meV, fio = 50meV, 
W = 20meV, Eo = 2eV. 

which is used in the following numerical calculations. The 
parameter A measures the strength of electron-phonon 
coupling, Qq is the mean phonon frequency, W defines 
the band width. Note that this correlation function takes 
into account multiple phonon emissions. A purely per- 
turbative treatment corresponds to the expansion of the 
exponential in Eq. @ to quadratic order in the coupling 
constant A. 

Fig. 121 shows the differential tunneling conductance 
dl/dV for different widths of the optical phonon band. 
The phonon satellites are most clearly visible in the case 
of completely localized Einstein phonons {W = ImeV). 
Each satellite then reproduces exactly the singularity of 
the quasiparticle peak in the superconducting DOS. The 
case of Einstein phonons seems to be excluded by the 
experimental data which shows broadened phonon satel- 
lites. Comparing qualitatively peak heights and peak 
widths of theory and experiment, we are led to the con- 
clusion that the phonon band has a width of roughly 
20meV. This width is probably due to the vibrational 
coupling of the apical oxygen atom to its neighbors. Since 
the electron-phonon coupling constant is quite large one 
might also expect broadening due to particle-hole excita- 
tions in the superconducting plane or due to the anhar- 
monicity of the crystal leading to a finite phonon lifetime. 

In Fig. 131 we show the dependence of the differential 
conductance on the electron-phonon coupling strength 
A. The stronger the coupling the more phonon satellites 
become visible in the spectrum. The amplitude of the 
nth satellite is smaller than that of the n — 1th satel- 
lite by roughly a factor of A^/n. A comparison of first 
and second satellite therefore allows for an estimate of 
the coupling strength A. From the experimental data we 
estimate A ~ 0.4. 

Fig. 21 evaluates the isotope effect in the second deriva- 
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FIG. 4: Isotope shift in the dPl/dV"^. Parameters are A = 
38meV, Q.o = 50meV, W = 20meV, A = 0.4, Eo = 2eV. The 
sign in the d^I /dV^ for negative voltages has been flipped. 

five (Pl/dV"^ of the current-voltage spectrum as mea- 
sured in the experiment 0. The isotope effect is most 
clearly visible for the first satellite which stems from 
a single convolution of the quasiparticle peak with the 
phonon spectrum. The second satellite stems from a 
twofold convolution with the phonon spectrum. Its sec- 
ond derivative d^I/dV^ is therefore much smoother, and 
the isotope effect less well seen. 

Experimentally , the lowest phonon frequency is ob- 
tained as the energy difference between the quasiparti- 
cle peak in the dl/dV (in our case at 38meV) and the 
maximum of the second derivative d'^I/dV^ (in our case 
at 78meV). This procedure is standard in conventional 
superconductivity where the strong singularity in the s- 
wave density of states guarantees a pronounced feature 
in the second derivative d'^I /dV^ . In the case of d-wave 
superconductivity this singularity is much weaker and 
consequently the maxima in the second derivative are 
less well defined. If the quasiparticle peak is weakened 
or if the second phonon satellite is investigated another 
procedure to determine the phonon frequencies can be 
preferable: Determining the positions of both maximum 
and minimum in the second derivative and taking their 
average one can directly estimate the central frequency 
of the phonon band. 

Inelastic tunneling processes show up as a pronounced 
peak in the d^I/dV"^ spectrum displaced by a phonon en- 
ergy from the main peak in the dl/dV associated with 
the DOS singularity which accompanies d-wave super- 
conductivity. The strength and form of the peak reflects 
those of the DOS singularity. This feature agrees with 
experiments of Lee et al. 0. A weaker signal also appears 
displaced by twice the phonon energies in agreement with 
experiment. The phonon energy, 50meV, chosen to agree 
with the experimental spectra, is close to the value ex- 
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pected for an apical 0-phonon and rather higher than the 
energy of the Big half-buckUng mode One point of 
divergence between the theory and the spectra reported 
by Lee et al. is the presence of a sharp minimum in the 
calculated (Pl/dV'^ spectra just above the peak. 

The STM experiments of Lee et al. show substantial 
spatial variations in the spectra with a clear correlation 
between the frequency shift of the side peaks and the 
local value of the superconducting energy gap. Such a 
correlation follows if the local bosonic mode couples di- 
rectly to the planar holes and if this electron-boson cou- 
pling is involved in the superconductivity. In our model, 
unlike those of Refs. |9l ll(i IllL Il2l|. such an explanation 
does not apply. However it is known empirically that 
the positions of the apical 0-ions influences the transi- 
tion temperature Tc of BSSCO superconductors [2^ and 
more generally that materials with longer Cu-apical O 
bonds have higher values of [13 ■ Therefore the anti- 
correlation between the local magnitude of the supercon- 
ducting gap and the local size of the frequency shift may 
reflect this empirical trend and is not incompatible with 
our model which ascribes these side peaks to inelastic 
tunneling processes in the barrier. 

Lcc et al. report an isotope shift when all 0-sitcs 
are occupied by O^^ vs. O^^. A selective isotope re- 
placement only on the apical 0-sites would discriminate 
between our proposal that the structure is due to apical 
0-phonons rather than due to the half-buckling mode 
involving displacements of planar 0-ions. 

In summary we conclude that inelastic tunneling pro- 
cesses involving the apical 0-ions are a viable explanation 
of the phononic structure in the STM tunneling spectra 
recently observed by Lee et al. and we propose selective 
isotope substitution as a test of our proposal. 
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